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One-dimensional delayed-detonation models of Type la supernovae: 
Confrontation to observations at bolometric maximum 
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ABSTRACT 

The delayed-detonation explosion mechanism applied to a Chandrasekhar-mass white dwarf 
offers a very attractive model to explain the inferred characteristics of Type la supernovae 
(SNe la). The resulting ejecta are chemically stratified, have the same mass and roughly the 
same asymptotic kinetic energy, but exhibit a range in 56 Ni mass. We investigate the con- 
temporaneous photometric and spectroscopic properties of a sequence of delayed-detonation 
models, characterized by 56 Ni masses between 0.18 and 0.81 M Q . Starting at 1 d after explo- 
sion, we perform the full non-LTE, time-dependent radiative transfer with the code CMFGEN, 
with an accurate treatment of line blanketing, and compare our results to SNe la at bolometric 
maximum. Despite the ID treatment, our approach delivers an excellent agreement to obser- 
vations. We recover the range of SN la luminosities, colours, and spectral characteristics from 
the near-UV to 1 fjm, for standard as well as low-luminosity 91bg-like SNe la. Our models 
predict an increase in rise time to peak with increasing 56 Ni mass, from ~ 15 to ~ 21 d, yield 
peak bolometric luminosities that match Arnett's rule to within 10 % and reproduce the much 
smaller scatter in near-IR magnitudes compared to the optical. We reproduce the morphology 
of individual spectral features, the stiff dependence of the 7£(Si) spectroscopic ratio on 56 Ni 
mass, and the onset of blanketing from Ti Il/Sc II in low-luminosity SNe la with a 56 Ni mass 
<0.3Mq. We find that ionization effects, which often dominate over abundance variations, 
can produce high-velocity features in Ca II lines, even in ID. Distinguishing between differ- 
ent SN la explosion mechanisms is a considerable challenge but the results presented here 
provide additional support to the viability of the delayed-detonation model. 

Key words: radiative transfer - supernovae: general 



1 INTRODUCTION 

All models for Type la supernovae (SNe la) invol ve the thermonu- 
clear disruption of a C/O white dwarf (WD) star dHovle & Fowler! 
Il960h . Fusion of 12 C and 16 to iron-peak elements provides the 
energy to unbind the star and accelerate the ejecta to expansion ve- 
locities on the order of 10000 km s -1 , while sy nthesizing < 1 Mq 
of 56 Ni to power the subsequent light curve dColgate & McKeg 
fl96^) . 

Nuclear combustion of C and O requires a large temperature 
(~ 10 9 K) in the progenitor star. At present, these conditions seem 
to be met when the WD mass approaches the Chandrasekhar mass, 
Mch ~ 1.4M0, through accretion of material from a compan- 
ion st ar (the so-called single-degenerate scenario; IWhelan & Iber] 
Il973l) . Alternatively, combustion may be initiated when two C/O 
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WDs coalesce (the d ouble-degenerate scenario; llben & Tutukovl 
ll984IWebbin3ll984l) . in which case the total mass of the sys- 
tem can differ significantly from Mch- While the single-degenerate 
scenario provides a natural explanation for the apparent homogene- 
ity of the SN la class, recent simulations have show n that violent 
WD m ergers are co mpatible with both sub-luminous jPakmor et alj 
l201Ch and normal jPakmor et aT]|2012l) SNe la. It is thus not clear 
today whether the Chandrasekhar mass represents a fundamental 
quantity for SNe la. 

In addition to multiple progenitor channels, several explosion 
mechanisms have been proposed. The propagation of the burn- 
ing front can proceed either sub-sonically (deflagration), or super- 
sonic ally (detonation). Pure detonations of C/O WDs (e.g.. lArriettl 
1 19691) bum the entire star to nuclear statistical equilibrium and fail 
to synthesize the intermediate-mass elements (IMEs ) at high ve- 
lociti es needed to reproduce SN la spectra (see, e.g., 



Branch et alj 



Pure deflagrations (e.g., iNomoto. Sugimoto & Ned 1 197 
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yield copious amounts of IMEs but correspond to less ener- 
getic explosions, with 56 Ni masses typically < O.4M0(see, e.g., 
iTravaglio et"ai]|2004h . and are thus considered a feasible mecha- 
nism only for the least luminous events. More importantly, signifi- 
cant amounts of unburnt C/O are expected to be mixed inwards to 
the WD centre (see, e.g.. lGamezo et al ] l2003l) . leading to the firm 
prediction of conspicuous lines of O I and C I at nebular times that 
have never been observed iKozma et all2005T) . 

To overcom e the shortcomings of both explosion mechanisms, 
iKhokhlovl ( fl99lh proposed the delayed-detonation model in which 
the explosion of a Chandrasekhar-mass C/O WD begins as a sub- 
sonic deflagration near its centre and then switches to a super- 
sonic detonation wave at some density, p tr , via a deflagration- 
to-detonation transition (or DDT). The initial deflagration pre- 
expands the WD and creates low-density conditions required for 
the production of intermediate mass elements. The ensuing detona- 
tion incinerates the remaining unburnt matter. By tuning p tr and the 
deflagration speed, one can synthesize varying amounts of 56 Ni and 
hence account for the observed variation in SN la peak luminosity, 
as well as reproduce the correct stratification of chemical elements 
in the ejecta inferred from spectroscopic observations. Additional 
variations in peak luminosity can be induced by va rying the C/O 
ratio and central density of the progenitor WD star l IHoflich et alj 

l2oToh . 

The delayed-detonation model has been success- 

fully tested against SN la observ ations in its ID (e.g., 

Hoflich. Khokhlov & Wheeled 1 19951) and multi-D (e.j 



Kasen. Ropke & Wooslevl l2009h versions. It also appears to be 



consistent with observations o f supernova remnants dBadenes et alj 
l200dl2008l : lFesen et alj|2007l). Using a g rid of 130 parameterized 



one-dimensional models, IWooslev et alj j2007l) studied the ejecta 
properties necessary to reproduce the fundamental properties of 
SN la light curves, and in parti cular the width-luminosity relation 
(Pskovskii 1977; Phillips 1993). They found a good agreement for 
ejecta with a total burnt mass of ~ 1.1 Mq (leading to comparable 
asymptotic kinetic energies in such models), provided the inner 
layers were strongly mixed, and highlight the adequacy of the 
delayed-detonation scenario in fulfilling these requirements. 
However, despite a significant effort in three-dimensional simula- 
tions (e.g., iGamezo. Khokhlov & Oranll2005l ; iRopke & Niemevej 
120071) . it is still unclear what may cause strong variations in 
transition density in the DDT mechanism. Moreover, the lack of 
unburnt carbon in delayed-detonation models is in conflict with 
observations of C II 658 A absorption in ~ 30% of SNe la with 
pre-maximum spec tra dParrent et alj 1201 it jThomas et alj l201ll ; 
iFolatelli et alj|2012l ; lsiiverman & Filippenkdl2012l) ~ 

Differentiating between these various models requires detailed 
radiative transfer simulations to match the observed evolution of 
SN la light curves and spectra form early to late times. At present, 
the observed properties of normal SN e la can be repro duced with 
detonations in single sub- Mr-h WDs (|sim et al.l l 2010l) and super- 
Mch double-WD mergers l lRopke et alJbOlaj Tin addition to the 
more conventional Chandrasekhar-mass delayed-detonation model. 
This may mean that observational differences between these mod- 
els are too small, or that uncertainties in the modeling approach are 
too large. 

Radiative-transfer modeling of SN la ejecta is complicated 
by the dominance of line op acity l IHoflich. Miiller & Khokhlov! 
ll993t|Pinto & EastmarfeOOObl) . with a mixed absorptive/scattering 
character that is often enforced, particularly in Monte-Carlo 
codes (see, e.g. Mazzali & Lucvlll993l ; iKasen. Thomas & Nugentl 
l2006t Isiml 12007^ " rather than computed directly as part of a 



non-Local-Thermodynamic-Equilibrium (non-LTE) solution. 
Other approaches treat some spe cies in full non-LTE, while oth- 
ers a re treated in LTE ( e.g., IHoflich. Wheeler & Thielemann! 
1 19981 : iBaron et all l2012h . Neglecting time-dependent ef- 
fects in the ra diation transport, which a r e known to be 
import a nt (see IHoflich & Khokhlovl 1 19961: I Pinto & Eastmar] 



| 2000bl: iDessart & Hillierl l20ld : Ijack. Hauschildt & BarorJ l2009j 
lHillier&Dessartl2012l) . numerous studies at bolometric maximum 
(and sometimes even later) impose a diffusive inner boundary, even 
though some spectral regions are no longer optically thick. This 
causes a systematic flux excess in the red, i.e., photons erroneously 
injected at the inner diffusive b oundary fail to fher malize and 
directly escape the ejecta (e.g., iMazzali et alj 120081) . The low 
densities and strong radiation field in SNe la is expected to drive 
level popu lations out of LTE , indicating the need for a full non-LTE 
treatment l lBaron et alJl 19961 ; Dessart et al., in prep.). 

As par t of a new and independent effort, we apply th e CM- 
FGEN code faillier&MiIieilll998l ; iHillier & Dessarlll2012h to the 
Chandrasekhar-mass delayed-detonation model of SNe la. CMF- 
GEN is a ID, non-LTE, time-dependent, radiative-transfer tool that 
allows for explicit treatmen t of non-thermal proces ses and non- 
local energy deposition (see IHillier & Dessart 2012 for a full de- 
scription of the code in the context of SN calculations). In this pa- 
per we focus on the photometric and spectroscopic properties of 
delayed-detonation models at bolometric maximum, to study how 
well this explosion mechanism can predict the observed diversity 
of SNe la at that well-defined time. 

In section[2]we present our grid of delayed-detonation models 
and radiative transfer modeling. In Sect. [3] we present our over- 
all strategy for comparing synthetic light curves and spectra. We 
then discuss the photometric (Sect. |4) and spectroscopic (Sect. [5} 
properties at bolometric maximum. We discuss comparisons to in- 
dividual SN la observations in Sect. [6] Discussion and conclusions 
follow in Sect. [7] 



2 NUMERICAL SETUP 

2.1 Delayed-detonation models 

A grid of one-dimensional delayed-deton ation mode l s was calcu- 
lated following the method described in IKhokhlovl dl99ll) . The 
starting point for all models is a Chandrasekhar-mass WD star 
(1.41 Mq) with a central density of 2.6 x 10 9 gem -3 , composed 
of equal mass fractions of 12 C and 16 (X ~ 0.5) with traces of 
22 Ne ( X 1.4 x 10~ 2 ) and so lar abundances for all other el- 
ements jAnders & Grevessell 19891) . The explosion hydrodynamics 
is described by the reactive flow Euler equations of fluid dynam- 
ics which are solved with a one-dimensional Lagrangian hydro- 
dynamics code. During the explosion, nuclear energy release and 
neutronization of matter were described using a five-equation re- 
action scheme for 12 C concentration, electron and mean ion mole 
fractions, nuclear binding energy, and the nuclear statistic equi- 
librium (NSE) progress variable. The reaction scheme takes into 
account the energy release during the major stages of explosive 
carbon-oxygen burning: (1) burning of carbon to oxygen, neon, 
magnesium and silicon, (2) burning of oxygen and establishment 
of nuclear quasi-equilibrium in Si-group nuclei, and finally (3) the 
onset of NSE during which the Si-group elements are transformed 
into iron-peak nuclei. The scheme takes into account neutroniza- 
tion and neutrino energy losses due to electron captures on iron- 
group nuclei. The neu tronization rates were updated according to 
ISeitenzahl etai] j2009h . 
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Table 1. Delayed-detonation model parameters and nucleosynthetic yields for selected species. 



Model p tr -Ekin f( 56 Ni) A/( 56 Ni) M(Ni) M(Co) M (Fe) M(Ti) M(Ca) M(S) M(Si) M(Mg) M(Na) Af(O) M(C) 

[gem- 3 ] [B] [kms- 1 ] [M Q ] [M Q ] [M ] [M Q ] [M ] [M ] [M ] [M ] [M ] [M Q ] [M Q ] [M Q ] 



DDCO 


3.5(7) 


1.564 


1.31(4) 


0.805 


0.857 


1.15(-2) 


0.160 


DDC6 


2.7(7) 


1.503 


1.21(4) 


0.727 


0.768 


9.10(-3) 


0.142 


DDC10 


2.3(7) 


1.485 


1.16(4) 


0.650 


0.689 


8.36(-3) 


0.137 


DDC15 


1.8(7) 


1.468 


1.12(4) 


0.558 


0.596 


7.87(-3) 


0.135 


DDC17 


1.6(7) 


1.410 


1.09(4) 


0.458 


0.495 


7.24(-3) 


0.131 


DDC20 


1.3(7) 


1.382 


1.04(4) 


0.344 


0.380 


6.43(-3) 


0.126 


DDC22 


1.1(7) 


1.301 


9.85(3) 


0.236 


0.271 


5.52(-3) 


0.120 


DDC25 


8.0(6) 


1.177 


8.36(3) 


0.180 


0.215 


4.37(-3) 


0.111 



1.19(-4) 2.62(-2) 0.101 0.160 3.51(-3) 5.88(-6) 5.29(-2) 1.18(-3) 

1.21(-4) 3.06(-2) 0.125 0.201 7.18(-3) 9.66(-6) 8.42(-2) 1.76(-3) 

1.24(-4) 3.68(-2) 0.149 0.234 9.55(-3) 1.20(-5) 9.67(-2) 2.11(-3) 

1.25(-4) 4.17(-2) 0.177 0.277 1.14(-2) 1.63(-5) 0.107 2.69(-3) 

1.25(-4) 4.45(-2) 0.200 0.314 1.65(-2) 2.43(-5) 0.139 3.77(-3) 

1.25(-4) 4.53(-2) 0.233 0.374 1.95(-2) 3.40(-5) 0.157 5.05(-3) 

1.24(-4) 4.04(-2) 0.255 0.432 2.45(-2) 6.07(-5) 0.186 8.14(-3) 

1.25(-4) 2.34(-2) 0.230 0.458 3.33(-2) 1.46(-4) 0.248 2.07(-2) 



Notes: Numbers in parenthesis correspond to powers of ten. The ratio of the deflagration velocity to the local sound speed ahead of the flame is a = 0.03 for 
all models; ptr is the transition density at which the deflagration is artificially turned into a detonation; Eki n is the asymptotic kinetic energy (units: 
IB = lBcthc = 10 51 erg); i>( 56 Ni) is the velocity of the ejecta shell that bounds 99% of the total 56 Ni mass (see Sect. 12.1 I f or a discussion). 



The equation of state includes contributions from electrons 
and positrons with an arbitrary degree of degeneracy and relativism, 
ideal ions, and equilibrium radiation. Screening of the 12 C+ 12 C 
nuclear reaction is taken into accoun t but Coulomb corrections t o 
the equation of state are ignored (see lBravo & Garcfa-Sen 
They change the relation between the initial central density and the 
WD mass by ~2-3% but do not noticeably affect the explosion 
predictions. The approximate kinetic scheme neglects differences 
between the quasi-equilibrium energy release in Si-group elements 
and the NSE energy release in Fe-peak elements and may overes- 
timate the energy generation rate during oxygen burning to silicon 
by up to ~20% as compared to detailed nuclear kinetics. 

The initial WD is assumed to be in hydrostatic equilibrium, 
and the discretized equations of energy and nuclear kinetics are 
solved simultaneously with high precision to prevent excitation of 
artificial pulsations during the WD explosion. Weak interaction 
rates a nd neutrino energy losses are taken from ISeitenzahl et alj 
j2009h . Detailed nucleosynthesis is calculated by post processing 
temperature, density, and neutronization histories of fluid elements 
with a detailed reaction network consisting of 144 individual nu- 
clei, ranging from protons and neutrons to zinc. Forward reaction 
rates, partition functions, and binding energies were taken from 
iFowler. Wooslev & Engeibrech.il Jl978h . lWooslev et all dl978b . and 
F. Thielemann (1993, private communication). Reverse reactions 
were calculated from the principle of detailed balance. 

The deflagration is propagated numerically with a prescribed 
velocity Ddcf = Qa s relative to unburnt matter, where a < 1 is 
a constant and a 3 is the local sound speed ahead of the deflagra- 
tion front. We use a — 0.03 in the models considered in this pa- 
per. The transition to a detonation is artificially triggered by rising 
Ddei above the maximum Chapman- Jouguet speed for a steady- 
state deflagration when the density ahead of the deflagration front 
reaches a prescribed transition density, p tr - This causes the onset 
of a self-sustained supersonic detonation wave after the accelerated 
deflagration front propagates ~ 3-5 computational cells. The cal- 
culations are run until 30-60 s past explosion, by which time the 
ejecta have expanded by two orders of magnitude and are in ho- 
mologous expansion. 

The earlier the transition from a deflagration to a detona- 
tion occurs (i.e. the higher the transition density ptr), the more 
efficiently the burning proceeds to large velocities with more 
56 Ni being synthesized at the expense of intermediate-mass ele- 
ments (Fig.Q}. From model DDCO to DDC25, the velocity of the 
ejecta shell that bounds 99% of the total 56 Ni mass drops from 



- 13000 kms -1 to ~ 8000 kms -1 , while the mass of synthe- 
sized 56 Ni goes from ~O.81M to ~0.18M Q (see Table Q). At 
the largest velocities one recovers the original composition of the 
progenitor WD star. 

Since the WD is almost entirely consumed in ID delayed- 
detonation models (there is at most ~ 0.02 Mq of unburnt carbon 
for model DDC25), the fraction of primordial C/O not burnt to 56 Ni 
is burnt to IMEs, with a comparable nuclear energy release in both 
cases. This leads to a unique property of the delayed-detonation 
model, whereby variations in the transition density results in SN la 
explosions with the same ejecta mass and roughly the same asymp- 
totic kinetic energy (Skin), but with varying amounts of 56 Ni. In 
our model set, -Ekin varies by only 50% for a change in 56 Ni mass 
by a factor of 4.5. 



2.2 Non-LTE time-dependent radiative transfer 

Our approach to solve the non-LTE time-dependent radiative- 
transfer problem for SN ejecta is versatile and generic. It ap- 
plies irrespective of the SN type and we thus use here the 
same approach as , e.g., in the recent study of pair-instability 
SN explosions by iDessart et al.1 d2012bl). The technique i s pre- 
sented in a sequence of papers dPessart & Hillierl 120051. [2008; 
ILL Hillier & Dessardl2o"l2l ; IDessart et alj|2012<j) . A thorough dis- 
cussion of the ap plication of the code CMF GEN to supernovae has 
been presented bv lHillier & Dessartl j2012h . 

Starting at 1 d after explosion, we simultaneously simulate the 
evolution of the full ejecta, including the gas thermodynamics state 
(non-LTE ionization and excitation, and temperature) and the radi- 
ation (either trapped or escaping). This evolution is controlled by 
cooling and heating processes. At early times, cooling stems pri- 
marily from expansion, which in SNe la corresponds to a factor of 
10 6-7 increase in radius from the original exploding white dwarf to 
the time of bolometric maximum at ~ 20 d. Heating is dominated in 
our simulations by the radioactive decay of 56 Ni and 56 Co nuclei. 

We solve the radiative-transfer problem using a non-LTE treat- 
ment. Hence, all atom/ion level populations are determined through 
a solution of the statistical equilibrium equations, coupled to the 
gas energy equations and the th and 1 st moments of the radiative- 
transfer equation. Level populations are ther efore controlled by ra- 
diative and collisional rates explicitly (see iHillier & Miller] 1 19981 
for details). We do not use the questionable assumption of LTE 
for the thermodynamical state of the gas, nor the Equivalent-Two- 
Level-Atom approximation for solving level populations. Together 
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Figure 1. Comparison of the mass fractions of Ni (top; 56 Ni shown as 
dotted lines), Si (middle), and O (bottom) in velocity space for our model 
set. The drop in 56 Ni mass fraction (e.g., at ~5000krns~ 1 in model 
DDC25) corresponds to the ejecta location where the deflagration transi- 
tions to a detonation. Note the absence of 56 Ni in the inner ejecta layers (up 
to v ~ 2000 km s _1 ) due to rapid electron captures. In the middle panel we 
overplot the location where the inward-integrated electron scattering opti- 
cal depth T es = 2/3 (dots), as well as minus the velocity at maximum 
absorption in the Si II 6355 A line (negative by convention; diamonds). 



with charge conservation, this yields a non-LTE ionization. Further- 
more, we allow for time dependence both for the radiative transfer, 
in order to compute accurately the transport of radiation through 
the ejecta, as well as in the statistical-equilibrium and energy equa- 
tions, in order to account for the time-depe ndent effects on the non- 
LTE io nization/excitation state o f the gas l lDessart & HillieJl2008l : 
see also lUtrobin & Chugaill2005l) . 

Opacity plays a critical role in SNe la. Through its impact 
on the ejecta optical depth, it controls in part the light-curve 



morphology dArnettl 1 19821 ; iHbflich. Miiller & Khokhlovl 1 19931 ; 
IPinto & Eastmanl2000bh . In the spectrum formation region, it mod- 
ulates the magnitud e of line blanketing ( IDessart 
iDessart et alj|2012bh . and is thought in S Ne la to lead to a sec- 
ondary maximum in near-IR light curves (Kasen 200^). In CMF- 
GEN, we incorporate line blanketing by explicitly treating bound- 
bound transitions between model atom/ion levels (this is part of the 
non-LTE solution). Hence, we use no ad-hoc prescription for the 
thermalization character of the gas, as typically done in the com- 
munity (sometimes in combination with LTE for the state of the 
gas). In the present simulations, we adopt a set of model atoms 
for which we obtain converged results (i.e., increasing further the 
size of the model atoms alters negligibly the resulting radiation and 
gas properties) — this required extensive testing, whose results will 
be presented in a separate paper (Dessart et al., in prep). In prac- 
tice, we include C I-IV, O I-IV, Ne I— III, Na I, Mg II— III, Al II- 

III, Si II-IV, S II-IV, Ar I— III, Ca II-IV, Sc II— III, Ti II— III, Cr II- 

IV, Mn II— III, Fe I-VII, Co II-VII, and Ni I I-VII, The number of 
levels (both super-levels and full levels; see iHillier & Millerll 19981 
for details) is given in the appendix, in Table Al. The highest ion- 
ization stages that are needed to model the transport at early times 
are no longer needed at bolometric maximum and thus these ions 
are omitted from that table. 

The output of the hydrodynamical calculations at 30-60 s past 
explosion is evolved to 1 d using a separate program. This program 
solves the energy equation given by the first law of thermodynam- 
ics, assuming the material is radiation dominated (which holds ev- 
erywhere apart from the o6 Ni-deficient higher-density inner ejecta 
below ~ 2000 km s -1 ); we allow for 56 Ni decay but neglect the 
diffusion of the associated heat. Given the high ejecta densities dur- 
ing the first day after explosion, and our focus on maximum-light 
spectra ~ 20 d later, these approximations are adequate. 

We then remap the resulting ejecta at 1 d into CMFGEN. The 
ejecta, whose velocity ranges from 500 to 45000 km s -1 and den- 
sity from~3x 10~ 9 to~3x 10 _ 15 g cm~ 3 , are typically resolved 
with 100-120 radial points, distributed evenly in optical-depth scale 
to have a minimum of 5-7 points per decade. For the initial re- 
laxation from the hydrodynamical input (see below), we keep the 
temp erature fixed and solve f or the full relativistic transfer in non- 
LTE dHillier & Dessar3l2012l) . The inner ejecta is very hot at that 
time (typically 10 5 K in the innermost 56 Ni-rich region), while the 
outer ejecta is at < 10 3 K. Since we model the full ejecta, this re- 
quires the simultaneous treatment of highly ionized species (which 
we limit to six-times ionized Fe/Co/Ni and five-times ionized O, 
for example) and neutral species with C I, O I, or Fe I. To avoid 
numerical overflow in this initial relaxation at 1 d (see details in 
IDessart & HillietolOl) . we have to enforce a floor temperature in 
the outer ejecta of 6000 K — this floor value is reduced to 2000- 
3000 K in subsequent time steps when we let the temperature vary 
in the time-dependent calculation. Together with the other approx- 
imations used for the evolution of the gas from < 100 s to 1 d, this 
implies that our simulations need a few time steps to relax from 
their initial conditions. 

We adopt a time increment of 10% of the current time un- 
til ~ 21 d past explosion, corresponding to 33 time steps from the 
initial time at 1 d. Subsequently, we use a fixed 2 d increment to en- 
sure a good time resolution. In the models with a large 56 Ni mass, 
the outer ejecta relaxes to the outward diffusion of heat released at 
greater depths, as well as radiative and expansion cooling, within 
a few 0. 1 d because of the close proximity of the 56 Ni to the outer 
ejecta. In the models with a low 56 Ni mass, this relaxation is slower, 
and takes ~ 1-2 d. None of these issues are important for the com- 
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parison to observations at bolometric maximum, ~ 15-21 d after 
explosion. Probably more relevant than any shortcoming of this ini- 
tial setup is the neglect of multi-dimensional effects, in particular 
associated with in ward mixing of unburnt C/O and o utward mixing 
of 56 Ni (see, e.g jGamezo. Khokhlov & Oranll2005h . 

We allow for non-local energy deposition from radioac- 
tive decay, using a Mo nte-Carlo approach for 7-ray transport 
dHillier & Dessartl |2012h . We also treat the non-thermal pro- 
cesses associated with the high-energy electrons produced by 
Compton scattering and photoelectric absorptio n of these 7-rays 
&Li, Hillier & Dessartll2012hlDessart et alj|2012al) . These processes 
are generally unimportant at early times and so we included them 
starting at ~ 15 d after explosion. 



3 STRATEGY 

In this study, we adopt a forward-modeling approach. We start 
with a grid of realistic explosion models and solve the time- 
dependent radiation-transport problem. An important validation of 
such delayed-detonation models is through direct confrontation to 
observations. We do not attempt to model a specific SN — in- 
stead we confront each of our models with a large database of 
SN la light curves and spectra. Our goal is to explore how well 
our range of delayed-detonation models (which span a factor of 4.5 
in 56 Ni mass) can reproduce the observed range of SN la proper- 
tie s, as for example de scribed by the different subclasses proposed 
bv lBranch etaTI J2006h . 

All models were evolved from 1 d after explosion until beyond 
bolometric maximum. Here we focus the discussion exclusively on 
the time of bolometric maximum, and defer to subsequent papers 
the study of the rise to peak, the early post-peak evolution and the 
width-luminosity relation, as well as the nebular phase. The spec- 
tra at each time step are integrated over frequency to yield bolo- 
metric light curves. We infer the time of peak bolometric lumi- 
nosity through a 3 rd -order polynomial fit to the light curve. The 
synthetic spectra are interpolated in time to yield maximum-light 
spectra for all models. Since we simultaneously compute spectra 
and light curves, we can directly compare SN luminosity and spec- 
troscopic morphology. For educational purposes, we can reveal the 
direct contributions of individual ions to the total flux by computing 
the synthetic spectra using the same converged radiative-transfer 
solution, but ignoring all bound-bound transitions of a given ion in 
the formal solution of the transfer equation (see Appendix [B}. 

In the next sections, we present the basic properties of our 
delayed-detonation model series and confront their photometric 
(Sect. [4]l and spectroscopic (Sect. [5} properties to the plethora of 
publicly-available SN la observations. Since the true bolometric 
luminosity is not known for most SNe la, we generate pseudo- 
bolometric (UVOIR) light curves from UB VRI magnitudes fo r 
both models and data using the method of IValenti et"aTl <2008l) . 
Observed magnitudes are corrected for extinction (in the MW as 
well as in the host galaxy) and placed on an absolute scale using 
the distance moduli inferred for each SN (Sect. [6). We then com- 
pare the synthetic and observed peak UVOIR luminosities, the ra- 
tio of which is noted Q uvo ir. We also compare the magnitudes and 
broadband colours at bolometric maximum in various optical and 
NIR bands. Observed spectra are de-redshifted and de-reddened 
and scaled to the absolute V-band magnitude derived from pho- 
tometry. To enable a better comparison of individual spectroscopic 
features, an additional scaling, F sc i, is applied to the data in order 
to match the mean synthetic flux in the range 4000-6000 A. In most 
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Figure 2. Ratio of the peak bolometric luminosity (ibol.peak) to the instan- 
taneous rate of decay energy (Z/d aaagj versus rise time to bolometric max- 
imum for our model set (see also lKhokhlov. Miiller & HoflicrJI 19931 their 
Fig. 29). Our results suggest that Arnett's rule (ibol.peak/^decay = 1) 
is accurate to within 10% for such delayed-detonation models. The inset 
shows the bolometric light curves up until a few days past bolometric max- 
imum (marked with a filled circle). 



cases Fsd ~ 1 to within 10%, and in all cases consistent with unity 
when taking into account errors on the adopted distances. The de- 
tailed comparison to individual observations achieving the "best- 
match" to each of our delayed-detonation models is presented in 
Sect. [6] 



4 PHOTOMETRIC PROPERTIES AT BOLOMETRIC 
MAXIMUM 

The fundamental asset of the delayed-detonation model is its ability 
to induce a variation in 56 Ni, controlled through the deflagration- 
to-detonation transition density. This introduces diversity in ejecta 
that have otherwise very similar credentials (mass, kinetic energy, 
energy lost to unbind the progenitor, energy lost to expansion). It is 
thus no surprise that this variation in o6 Ni drives the main contrast 
in model results. 

The first observable affected by the variation in Ni mass is 
the peak luminosity, which correlates with the heating rate from de- 
cay. In the sequence DDC25 to DDC0, the peak bolometric lumi- 
nosity ranges from 3.64x 10 42 erg s -1 to 1.77x 10 43 erg s" 1 (i.e., 
from 9.5x 10 8 L to 4.6x 10 9 L ). We obtain the following rela- 
tion for our model set between the peak bolometric luminosity and 



the 56 Ni mass: 



^bol,pcak 



2.21 x 10 43 x M(J 



; Ni) 



M. 



erg s 



The degeneracy of progenitor and explosion properties, varying pri- 
marily in the amount of heating through decay, is in part the cause 
of this linear relation to 56 Ni mass. 

For all models the ratio of the peak bolometric luminosity to 
the instantaneous rate of decay energy (noted Q 7 in Table |2]l is 



within ~ 1 0% of unity ( Fi 
nett rule" jArnettl \197 



, thus supporting the so-called "Ar- 
19821) . This rule proposes that the light 
curves of the bolometric luminosity and the instantaneous decay 
rate should cross around the time of bolometric maximum. At this 
time, the ejecta Rosseland-mean optical depth in our simulations 
is a few tens, so that there is still trapped radiation within the SN, 
making the optical light curve exceed the 7-ray light curve after 
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Table 2. Various quantities at bolometric maximum characterizing our delayed-detonation model set. 
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ri\/T.~,i 

UVIqJ 


LKm s 


[1U JS.J 
















DDCO 


0.805 


15.71 


1.77(43) 


0.91 


-0.17 


0.848 


9699 


1.527 


8.46(-l) 


1.18(-1) 


1.58(-5) 


5.58(-ll) 


2.52(-2) 


4.00(-2) 


2.76(-7) 


DDC6 


0.727 


16.89 


1.62(43) 


0.98 


-0.13 


0.809 


9297 


1.422 


7.24(-l) 


1.15(-1) 


2.14(-5) 


8.63(-ll) 


4.00(-2) 


8.95(-2) 


6.83(-7) 


DDC10 


0.650 


17.44 


1.45(43) 


1.01 


-0.09 


0.793 


9129 


1.342 


6.13(-1) 


1.08(-1) 


2.38(-5) 


1.10(-10) 


4.89(-2) 


1.36(-1) 


1.12(-6) 


DDC15 


0.558 


17.76 


1.22(43) 


1.03 


-0.02 


0.786 


9019 


1.263 


4.54(-l) 


9.27(-2) 


2.43(-5) 


1.48(-10) 


5.60(-2) 


2.14(-1) 


2.05(-6) 


DDC17 


0.458 


18.90 


1.03(43) 


1.10 


0.05 


0.738 


8647 


1.197 


3.68(-l) 


8.69(-2) 


2.43(-5) 


1.60(-10) 


5.94(-2) 


2.56(-l) 


2.52(-6) 


DDC20 


0.344 


19.38 


7.65(42) 


1.11 


0.15 


0.726 


8504 


1.091 


1.95(-1) 


6.70(-2) 


2.02(-5) 


2.10(-10) 


5.65(-2) 


3.62(-l) 


4.46(-6) 


DDC22 


0.236 


19.52 


5.10(42) 


1.09 


0.65 


0.721 


8306 


0.980 


8.71(-2) 


5.09(-2) 


1.50(-5) 


2.69(-10) 


4.68(-2) 


4.61(-1) 


7.69(-6) 


DDC25 


0.180 


20.50 


3.64(42) 


1.06 


0.96 


0.626 


7338 


0.920 


3.49(-2) 


4.06(-2) 


9.58(-6) 


3.44(-10) 


3.29(-2) 


5.51(-1) 


2.18(-5) 



Notes: Numbers in parentheses correspond to powers of ten. Q 7 is the ratio of the peak bolometric luminosity to the instantaneous rate of decay energy. The 
Lagrangian mass coordinate M T , velocity v T , temperature T T , and mass fractions X are given at the location where the inward-integrated electron-scattering 
optical depth r cs = 2/3. Note that the reported 56 Ni mass fraction corresponds to a time immediately after explosive burning ceases, rather than the time of 
bolometric maximum. 



peak. It is only at nebular times that the SN radiation is directly tied 
to the instantaneous energy deposition rate. A fraction of the decay 
energy will inevitably escape the ejecta, even as early as maximum 
light. In our model series, the energy actually deposited in the ejecta 
is 85%-95% of the dec ay energy, resulting in a 4% -17% higher ef- 
fective Q 7 (noted Q in lHoflich & Khokhlov|[T99rl . 

We also o btain a fair agreement with t he delayed-detonation 
model results of lHoflich & Khokhlovl i ll 9961) for the peak bolomet- 
ric luminosities, but our rise times are 3-8 d longer for a given 
56 Ni mass and show a systematic increase for lower 56 Ni mass. 
The origin of the differences may be the opacity (and its associ- 
ated impact on optical depth and the diffusion time), whose influ - 
ence has been studied, e.g., bv lKhokhlov. Miiller & Hofiici] ( ll993l) : 
IPinto & Eastmanl l feOOOalfbT) . Since the pre-peak phase is not the fo- 
cus of this study, we defer a detailed discussion of this result, and 
comparisons to alternate works, to Blondin et al. (in prep). 

At bolometric maximum, the "photospheric" radijj are on the 
order of ~ 1.3— 1.4x 10 15 cm irrespective of the model, so that the 
range of luminosities reflects a modulation in temperature, from 
9200 K in model DDC25 to 15270 K in model DDCO at this loca- 
tion. Consequently, a trend of higher ejecta ionization is obtained 
for models with increased o6 Ni (Fig. [3}. Combined with abundance 
variations in the spectrum formation region, these variations are at 
the origin of the basic diversity of SN la maximum-light spectra 
(Sect.HJ. 

With our sample of delayed-detonation models, which spans 
a range in 56 Ni mass from 0.18 to 0.81 Mq, we reproduce the full 
range of observed luminosities, 73-band magnitudes and colours 
of SNe la at maximum light (Fig. |4j- For example, our model 
set covers B — R values from —0.17 to 0.96 mag at bolomet- 
ric maximum (Table |2). and misses only a few of the brightest 
and bluest SNe la at peak (most of these are 91T-like). We also 
reproduce the small luminosity scatter in the near-infrared bands 
with respe ct to the optical (see Table CI), as fir st found observa- 
tionally by|Kr tsctunas, Phillips & Suntzefdd2004) . The remarkable 
achievement of the delayed-detonation mechanism, which is not 
new, is that a single parameter can explain the most fundamental 



1 In SNe la, there is no clearly-defined photosphere due the ubiquitous 
presence of lines and the sub-dominance of absorption/continuum pro- 
cesses. Hence, by this term, we broadly refer to the spectrum-formation 
region, and adopt for that region the location where the inward-integrated 
electron-scattering optical depth is 2/3. 




O 
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M( 56 Ni) [M ] 

Figure 3. Ionization fraction at bolometric maximum for selected species as 
a function of 56 Ni mass, at the location where the inward-integrated elec- 
tron scattering optical depth r cs = 2/3. We define the ionization fraction as 
iX t+ 1 X tJr , where X l+ is the mass fraction of ionization stage i 
for element X. Numbers below the curve for the oxygen ionization fraction 
correspond to the DDC model number. 

properties of SN la light curves. In the next section, we show how 
well it reproduces the fundamental spectroscopic properties of a 
wide range of SNe la at bolometric maximum. 



5 SPECTROSCOPIC PROPERTIES AT BOLOMETRIC 
MAXIMUM 

5.1 Spectral Diversity, ionization sequence, and line 
blanketing 

The spectral diversity of these delayed-detonation models, asso- 
ciated with variations in ejecta characteristics (temperature, ion- 
ization) and photometric properties (luminosity, colour) reflects 
closely the observed diversity of SNe la at maximum light (Fig.[5j 
see also Sect.|6]l. This is a real achievement given the much greater 
complexity of spectra compared to photometry (for example, a 
given model may reproduce t he colour of a SN l a but may match 
the spectrum poorly; see, e.g.. lBlondin et al 120T l]). The theoretical 
spectral-energy distributions (SEDs) match satisfactorily the entire 
observed wavelength range from the blue, luminous, but spectro- 
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Figure 4. Intrinsic peak Mb (left) and UVOIR luminosity (right) vs. intrinsic B — R colour at B-band and UVOIR maximum, res pectively, for our seq uence 
of delayed-detonation models. We overplot (grey crosses) the observed values using a wide sample of SNe la (similar to that used in lBlondin et al.l201 lh. Star 
symbols correspond to objects that will be discussed individually and in greater detail in Sect. [6] Synt hetic a nd observed UVOIR luminosities were derived 
from UBVRI magnitudes (only BVRI in the case of SN 2007hj) using the method of lValenti et all l2008t) . We provide tabulated values of the bolometric 
luminosity and U BVRIJHK magnitudes at bolometric maximum for our models in Table C 1 . 



scopically normal SN 2009ig (model DDCO), to the redder, low- 
luminosity 91bg-like SN 1999by (model DDC25). For each in- 
termediate model, we also find suitable matches. Models with a 
larger 56 Ni mass achieve a higher peak bolometric luminosity and 
bluer colours, which in spectra is best reflected with the stronger 
flux shortward of ~ 4000 A. All these properties directly reflect the 
varying levels of temperature/ionization discussed in the preceding 
section. 

With a total 56 Ni mass of 0.81 M , model DDCO matches 
the peak luminosity of several 91T-like SNe la but fails to repro- 
duce their blue spectra characterized by lines of Fe III and very 
weak Si II 6355 A absorption. In the future, we will investigate if 
a delayed-detonation model with a (slightly) larger 56 Ni mass than 
DDCO would yield a better agreement with 91T-like objects. 

Irrespective of 56 Ni mass, our models are systematically faint 
in the UV. The total integrated flux blueward of 2500 A represents 
between -0.6% (model DDCO) and < 0.01% (models DDC22 
and DDC25) of the bolometric flux. This agrees with the few 
observations shown in Fig. [5] a subset of which have UV pho- 
tometry obtained w ith the Hubb le Space Telescop e (SN 1992 A, 
lKirshneretaLlll993l ; SN 2004dt IWang et alj|2012l) a nd the Swift 
satell ite (SN 2005cf, IWang et all l2009bl : SN 2009ig, iFolev et alj 
120121) . In addition, the inferred uvw2 — v colour jf] range between 
~ 2.5 mag (DDCO) and ~ 5 mag (DDC25), in line with those in- 
ferred by dMilneetalJl2010l) using Swift observations of a large 
sample of SNe la. 

The key spectroscopic signatures of SNe la are reproduced, 
such as the near-constant strength of the distinctive Si II 6355 A 
line in observations of low to high luminosity SNe la. In con- 
trast, the relative depth at maximum absorption of the Si II 5972 A 
line is a strongly varying function of peak luminosity and colour 
(see Sect. |5.2| l. We also predict the presence at maximum light of 
the Ca II near-infrared triplet (8498, 8542, and 8662 A) and the 
emergence of a strong O 1 7773 A line for SNe la with less than 
~ 0.4 Mq of 56 Ni. Both features are discussed in greater detail in 
Sectionsl531andl54l 



2 The uvw 2 and v filters have a central wavelength of 1941 A and 5441 A, 
respectively terown et al.l20ld) . 



While the colour of the SED is in part set by the tempera- 
ture/ionization of the spectrum formation region, another critical 
ingredient is line blanketing. At the peak of the light curve, the 
bulk of the radiation emerges from regions where the ejecta com- 
position is dominated by IMEs and IGEs (see Tableland Fig. [TJ, 
whose associated line opacity supersedes the contributions from 
electron scattering and continuum processes. These different opac- 
ity sources may overlap but often do not, so that they effectively 
blanket the entire spectrum, leaving no gaps for radiation to escape 
except at much longer wavelengths. Hence, blanketing induces sig- 
nificant redistribution of the flux to longer wavelengths through ab- 
sorption and fluorescence. 

Figure [6] illustrates the impact of lines on the maximum-light 
spectra and suggests two types of influence. First, Si II, Ca II, 
Mg II, or O I produce few strong and isolated lines through transi- 
tions between low lying states, which affect a limited wavelength 
range, although for Ca II H&K, the width of the absorption can be 
as broad as the U -band filter (model DDC20; see Fig. B6)! Second 
and more dramatic is the presence of forests of lines from some 
species/ions, which effectively blanket the SN radiation. For ex- 
ample, S II causes significant blanketing in the 5500 A region (the 
effect is largest for models DDC15-DDC17) — the associated lines 
overlap and prevent the identification of any of them individually 
(although the strongest transitions produce the characteristic "W" 
feature; see Appendix iBl. Ionization also tunes the effects of Fe 
opacity, which is primarily due to Fe III at high 56 Ni mass and to 
Fe II for models less luminous at peak than DDC15. 

More striking still is the sudden increase of Ti II and Sc II 
opacity in models with < 0.4 M of 56 Ni (DDC20-DDC25), due 
to a shift to a lower ionization state (Fig. [3j. This largely explains 
the progressive drop in the emergent flux blueward of ~ 4500 A in 
our synthetic spectra, which mirror the observations, as we progress 
in the model sequence from DDCO to DDC25, i.e., from blue lu- 
minous SNe la to red low-luminosity 91bg-like SNe la. This effect 
may seem surprising given the low abundance of Ti and Sc in the 
ejecta (Table[T](. The distribution of Ti is quite uniform throughout 
the ejecta with a mass fraction of ~ 10~ 5 at r cs = 2/3 and beyond, 
such that the increasing impact of Ti II is entirely attributed to an 
ionization effect. In addition, abundance variations play a role for 
the Sc II lines, since the Sc mass fraction is ~ 10~ 10 at r cs = 2/3 
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Wavelength [A] 

Figure 5. Comparison of synthetic spectra of delayed-detonation models (black) to observed SNe la near maximum light (red). The tickmarks on the ordinate 
give the zero-flux level. The time label corresponds to days from pseudo-bolometric (UVOIR) maximum. The observed spectra have been de-redshifted, 
de-reddened, and scaled to match the absolute V-band magnitude inferred from the corresponding SN la photometry. An additional scaling (F sc \; green label) 
has been applied to the observed spectra t o reproduce the me an synthetic flux in the range 4000-6000 A. For SN 2004dt we add the low-resolution HST ACS 
grism spectrum (dotted line) published by Wan g et al.l J20121) . 
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Figure 6. Comparison of synthetic spectra in our model series at bolometric maximum illustrating individual ion contributions to spectral features. Besides 
the total spectrum (black), the left panel shows the influence of O I, Mg II, Si II, S II, Ca II, and Sc II; the right panel shows that of Ti II, Co II, Co III, Fe II, 
and Fe III — the coloured lines correspond to synthetic spectra computed ignoring all bound-bound transitions of the corresponding ion. For better visibility 
of weak flux variations in the red, we show the quantity X 2 F\. 



(see Table O, but reaches ~ 10~° at larger velocities. The influ- 
ence of Ti Il/Sc II opacity becomes particularly severe starting at 
model DDC22, which marks a transition to 91bg-like spectra, char- 
acterized by enhanced flux blocking over the range 4000-4500 A. 
Note how this transition occurs over a very narrow range in 56 Ni 
mass (from 0.34 M for DDC20 to 0.24 M for DDC22), such 
that in principle a strict upper limit on the 56 Ni mass can be in- 
ferred from spectra with a 91bg-like morphology (in the context of 
a Chandrasekhar-mass delayed-detonation model). 



5.2 The Si II 6355 A line and the fc(Si) ratio 

The defining spectroscopic feature of SNe la, namely the 
blueshifted Si II 6355 A doublet absorption, remains broad 
(FWHM~ 12000 km s _1 ) in models DDC0-DDC17 and only 
becomes noticeably narrower from DDC20 onwards, with 
FWHM~8000kms _1 in model DDC25, reflecting the stronger 
confinement of Si to the inner ejecta in models with a lower 
deflagration-to-detonation transition density (Fig.Q3. The velocity 
at maximum absorption (w a bs) approximately reflects that of the 
peak Si abundance (Xsi ~ 0.6, regardless of the 56 Ni mass; see di- 
amonds in Fig. [l}, and varies smoothly between ~ — 17000 km s -1 
for DDCO and ~ -9000 km s" 1 for DDC25. 

The relative depth at maximum absorption varies little 
amongst our model series, as the line is optically thick at this lo- 



cation in all models (the Sobolev optical depth for the transition 
at 6347 A varies between ~ 2 for DDCO and ~ 387 for DDC25). 
This contrasts with the neighbouring Si II 5972 A line which be- 
comes progressively stronger with decreasing 56 Ni mass, transi- 
tioning from optically thin to thick (rs b(5979 A) > 1 at maxi- 
mum absorption for models DDC20-25) a nd reflecting the drop i n 
the Si ionization fraction (Fig. |3j see also lHachinger et alj|2008h . 
The relative behaviour of both lines is a well-known spectroscopic 
characteristic of SN e la, and quantified in terms of the 7?.(Si) ratio 
dNugentetaklll995h . In our models 7?.(Si) varies between < 0.1 
(DDCO) and ~ 0.8 (DDC25), in good agre ement with the observed 
variation (Fig. [7J. Contrary to claims by iGarnavich et alj J2004I) 
we find no evidence for Ti II absorption in this wavelength region 
that could bias the 7£(Si) ratio for low-luminosity SNe la (see Ap- 
pendix |B). 



iBranch et al. IHjoS) use the pseudo-equivalent widths (pEW) 
of the Si II 6355 A and 5972 A lines to classify SNe la into four 
subgroups, organized around the "Core Normal" subgroup which 
occupies a central position in this two-dimensional parameter 
space. SNe la with smaller and larger pEW(6355 A) are termed 
"Shallow Silicon" (includes 91T-like SNe la) and "Broad Line", 
respectively, while those with larger pEW(5972A) are dubbed 
"Cool" (includes 91bg-like SNe la), in reference to their lower 
ejecta temperatures. Our models belong to only two subclasses 
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model set. We ov erplot (grey crosses ) observed values from a large sam- 
ple of SNe la (see lBlondin et al.l2012h . Star symbols correspond to objects 
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of this scheme, namely Broad Line (models DDCO-17) and Cool 
(DDC20-25) SNe la. 

The classification scheme o f IWang et alj j2009ah is based on 
the velocity at maximum absorption of the Si II 6355 A feature, 
and divides SNe la into Normal and High- velocity subclasses. 
SNe la with 91T-like or 91bg-like spectra (including SNe la with 
7£(Si) > 0.5) are not classified according to this scheme. As 
noted in Sect. 15.1 1 there are no 91T-like spectra in our model se- 
ries, while models DDC20-25 either have TZ(Si) > 0.5 (DDC20) 
or display 91bg-like spectra (DDC22-25). The remaining models 
(DDCO-17) would all be part of the High-velocity su bclass (i.e., 
«abs < -12000 kms -1 ; see, e.g jBlondm et alj|2012h. 

There are thus no C ore Normal branch et alj 2006t) or Nor- 
mal l lWang et alj 12009 i) SNe la amongst our model set. The 
Si II 6355 A line in our synthetic spectra is systematically broader 
and more blueshifted than observed in these two subclasses, sug- 
gesting the ejecta kinetic energy in our models is too large to ac- 
count for the bulk of SNe l4j. However, the impact is often only 
visible in the Si II lines (4130 A, 5972 A, and 6355 A), while other 
spectroscopic features are well matched (e.g., Ca IlH&K and the 
S II "W" feature in the comparison between model DDC10 and 
SN 2005cf in Fig.|5]l, such that some of our models closely resem- 
ble Core Normal SNe la despite the Si II 6355 A mismatch (DDCO 
and DDC10). If this were solely a problem related to the ejecta 
kinetic energy then all lines should be too broad. Our limitation 
to ID also overlooks potential viewing-angle effects in an aspher- 
ical ejecta, which can significantly im pact the width and bl ueshift 
of the Si II 6355 A line, as shown by iBlondin et alj 1 201 1 ) based 
on the 2D delayed-detonation models of lKasen. Ropke & Wooslevl 
<2009h . 



5.3 The O 1 7773 A line as an abundance tracer 

The evolution in strength and morphology of the O 1 7773 A line 
profile in our synthetic spectra reproduces well the observed se- 
quence. This line becomes progressively stronger with decreasing 



3 The simplified reduced nuclear kinetics scheme used in our models (see 
Sect. |2. U may lead to an overestimation of the asymptotic kinetic energy. 
Hydrodynamical simulations with detailed nuclear kinetics are planned to 
quantify the difference. 



56 Ni mass along our model sequence (Fig. [6] left; green line). 
It is essentially non-existent in models DDCO- 10, barely notice- 
able in models DDC15 and DDC17, only to strongly imprint the 
spectra from model DDC20 onwards. This evolution follows the 
gradual decrease in the ionization fraction for O from ~ 2.0 to 
1.0 (Fig. [3}, but seems primarily influenced by the varying oxy- 
gen abundance with 56 Ni mass (factor ~ 5 variation in total oxy- 
gen mass between DDCO and DDC25) and the resulting chemi- 
cal stratification in velocity. The velocity at which X(0) = 10 -3 
drops from ~ 17500 kms -1 in model DDCO to ~ 9200 krns" 1 
in model DDC25, while the corresponding "photospheric" veloc- 
ities are 9699 and 7338 km s -1 , respectively (Table O. We thus 
associate this behavior of the O I line with an abundance effect re- 
sulting from the ejecta stratification in velocity. 

The O 1 7773 A line can in principle be used as a probe of 
left-over unburnt oxygen. However, unlike C II lines which pro- 
vide an unambiguous signature of unburnt carbon, the only definite 
signature of unburnt oxygen (as opposed to O produced via carbon 
burning) in our models is through the detection of absorption out 
to velocities where we recover the initial X(0) ~ 0.5 composi- 
tion. The minimum velocity at which X(0) ~ 0.5 varies between 
~ 20000 km s" 1 (model DDC25) and ~ 35000 km s~ x (model 
DDCO), and we obtain no O 1 7773 A absorption out to these veloc- 
ities in any of our models. This is likely a shortcoming of spherical 
symmetry, as multi-D effects are considered essential to preserve 
O (and C) unburnt in the ejecta, even at the relatively low veloc- 
ities correspon ding to the spectral-formation reg ion at bolometric 
maximum (see[G amezo, Khokhlov & Oranll2005() . 



5.4 Impact of Ca II lines and high-velocity features 

The large drop in flux around ~ 3500 A between model DDC17 
and DDC20 is mainly due to an opacity increase associated with 
Ca II. While the Ca ionization fraction at r os = 2/3 stays fairly 
constant from models DDCO to DDC25, it changes significantly 
at velocities > 15000 kms -1 (Fig.[8}. Since the Ca IlH&K lines 
remain optically thick out to large velocities, absorbing well be- 
yond 30000 km s _1 in all models, such ionization changes have a 
noticeable impact on the emergent flux in the U band. 

The effect is more clearly visible in the Ca II near-infrared 
triplet, which remains weak in all models up to DDC17 and be- 
comes progressively stronger up to DDC25, with a line-profile mor- 
phology displaying high-velocity components (see inset in Fig. [8}. 
Such high-velocity features ( HVFs) were associate d with density or 
abundance enhancements bv lMazzali eTai] ( 120051) . possibly result- 
ing from multi-dimensional effects in the explosion itself or from 
interaction with circumstellar material. Our simulations suggest 
that HVFs in Ca II can even arise in spherically symmetric ejecta, 
primarily through optical-depth effects due to ionization changes 
occurring at large velocity. Abundance effects play no role here: 
the Ca mass fraction at 15000 km 1 is three orders of magnitude 
lower in model DDC25 (~ 4.1 x 10~ 5 ), which displays promi- 
nent HVFs, compared to model DDCO (~ 3.4 x 10" 2 ), which dis- 
plays a single weak absorption. Furthermore, there is no density 
bump at these velocities. Obviously, when such HVFs are associ- 
ated with a va riation in the polar ization signal across the line (e.g., 
in SN 2001eh lWang et al]l2003l) . a departure from spherical sym- 
metry must occur, but these two signatures (HVF and polarization) 
can arise independently of one another. 

In all of our models with more than 0.4 M of 56 Ni (DDCO- 
17), the Ca II H&K absorption feature is strongly contaminated by 
the Si II 3858 A doublet, resulting in double-absorption features 
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Figure 8. Ionization fraction for Ca at bolometric maximum as a function 
of velocity. The solid dots indicate the location where the inward-integrated 
electron scattering optical depth t cs = 2/3. The inset shows Ca II near- IR 
triplet line profiles computed as the ratio of the full spectrum (.Fa, All) to 
that excluding all Ca II bound-bound transitions (Fa, leas), revealing high- 
velocity absorption features in models with less 56 Ni. 

(see Appendix [Bj. These are occasionally mistaken for "photo- 
spheric" a nd "high-velocity " components of Ca IlH&K (see dis- 
cussion in lLentz et alj2000f) . 



6 INDIVIDUAL COMPARISONS OF 

DELAYED-DETONATION MODELS TO 
OBSERVATIONS AT MAXIMUM LIGHT 

For each Chandrasekhar-mass delayed-detonation model in our set 
we are able to find a good match to an observed SN la (Fig.O. 
In most cases there are several observed SNe la that match a given 
model — this is not surprising given the overall s pectral homo- 
genei ty of SNe la with a given luminosity (see, e.g.. lBlondin et alj 
120121) . In those cases we preferentially select SNe la whose spec- 
tra cover the largest wavelength range (i.e. extending blueward of 
Ca IlH&K and redward of the Ca II NIR triplet). The compari- 
son between synthetic and observed spectra is largely based on a 
"chi-by-eye" estimate of the overall SED match, with particular 
emphasis on the flux level in the U and B bands, as well as in- 
dividual spectroscopic features. This approach is fine to assess the 
broad compatibility of delayed-detonation models — given the rec- 
ognized influence of multi-D effects and our ID approach, a perfect 
match may be coincidental. We use the Supernova Identification 
(SNID) code of lBlondin & Tonrvl ( 120071) to ease our exploration of 
the largest possible data set, but make no attempt at quantifying the 
strength o f the correlation betw een synthetic and observed spectra 
as done bv lBlondin et al. 

6.1 DDCO (M 56Ni = 0.81 Mq) vs. SN 2009ig 

The spectra of SN 2009ig frolev et alj|2012h are well reproduced 
by DDCO over the wavel ength range 2000- 10000 A. We use the 
same distance modulus as iFolev et alj J2012h (from lTullv & Fished 
Il988l) . but a larger extinction (Ay = 0.3 mag from our own 
fits to the light curves using the BayeSN stat i stical model of 
iMandel. Naravan & Kirshneil l201ll) . IFolev et alj fe012h do note 
however the presence of significant interstellar Na I absorption 



(EW = 0.4 A) at the redshift of the host galaxy, c onsistent with 
visua l extinctions in the range < Av < 1 mag (see lBlondin et alj 
l2009l their Fig. 5). 

SN 2009ig is a "Core Normal" SN la according to th e 
iBranch et al] J2006h classification scheme dBlondin et alj I2012T ). 
We are thus unable to reproduce its narrow Si II 6355 A line pro- 
file (see Sect. 15. 2\ . There is only a minor mismatch, however, in the 
width of other features, such as the Ca II H&K absorption (contam- 
inated by Si II 3858 A in the blue; see Sect.[5l]and Fig. Bl). The 
comparison at —3.3 d is impressive overall, although it is not as 
good for spectra closer to maximum light, in particular in the region 
3000-3500 A where the model is brighter than observations. The 
ratio of peak UVOIR luminosities is comparable (Quvoir = 1.12) 
to the scaling to match the mean fluxes in the 4000-6000 A region 
(-Fsci — 1-03), confirming the good correspondence between the 
model and the data. 

The peak lumi nosity and colours of DDCO match well th e lu- 
minous SN 1999ac dGaravini et alj|2005l ; IPhillips et alj|2006l) and 
SN 2001 eh l lSauer et alfeOOSt) . although our synthetic line profiles 
are systematically too broad (in particular Ca II H&K, in addition 
to Si II 6355 A) with respect to the observations. 

Model DDCO also reproduces the overall SED and spec- 
troscopic fe atures of the lumino us yet spectroscopically normal 
SN 1999ee faamuv et alj |2002|) . and matches well its peak lu- 
minosity (Quvoir = 0.98), although the flux level blueward of 
~ 3500 A is too high in our synthetic spectrum. 

6.2 DDC6 (Ms6 Ni = 0.73 M©) vs. SN 2004dt 

SN 2004dt has some interesting properties. Its Si II 6355 A line dis- 
plays pro minent HVFs up until a few days past maximum light, 
which led lAltavilla et alj |2007) to invoke an enhanced density 
structure in the outer layers. Furthermore, it displays a UV ex- 
cess around maxim um light compared to other SNe la, which led 
I Wang etalj J2012I) to invoke viewing-angle effects associated with 
large-scale ejecta asymmetries. Interestingly, our smooth ID model 
DDC6 reproduces the line widths, the overall SED, and the strong 
UV flux of SN 2004dt. (The dotted line in Fig. [5] shows an HST 
UV grism spectra taken with the ACS instrument — the resolution 
is very low so we are not concerned with small-scale absorption 
features.) 

We find an excellent match to the peak luminosi t y, col ours, 
and early-time spectra of SN 2006cp dBlondin et alj |2012|) . al- 
though the published spectra span a limited wavelength range 
(3500-7500 A) and none are within 5 d from B-band maximum 
light. 

6.3 DDC10 (A/ 5 6 Ni = 0.65 M Q ) vs. SN 2005cf 

SN 2005cf is a "go lden standard" amongst SNe la according to 
IWang et alj <2009bh . The match of the maximum-light spectra to 
model DDC10 is excellent overall in the range 2000-10000 A, al- 
though the Si II lines (4 130 A, 5972 A, and 6355 A) are signifi- 
cantly narrower and less blueshifted than in our synthetic spectrum 
(as expected since this SN is part of the Core Normal subclass; 
see Sect. 15.2b . The widths of other features, such as Ca II H&K 
and the S II "W", are nonetheless compatible with our model. The 
Ca II H&K absorption displays a double-absorption profile whose 
blue component is mostly due to contamination by Si II 3858 A 
(Sect. [51] and Fig. B3). Model DDC10 also reproduces the UV 
flux for this SN, measured using the Swift satellite, and matches 
perfectly the observed UVOIR flux (Q UV oir = 1.00). 
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Table 3. Comparison of delayed-detonation models to observed SNe la near maximum light. 



Model 


1V1 ^ IN 1 J 


•rise 


^bol 


r 

^uvoir 




Ker. 


A* 


ivietnoa 


rj\D — V ) 


uvoir 


Ai 

^-^uvoir 


-^scl 




TM^ 1 
L ivl J 


laayj 


[ergs 


[erg s 1 ] 






[mag] 




[mag] 




laayj 




(1) 


(Z) 


(3) 


W 


(?) 




(I) 


/ Q * 
\°) 


(V) 




U 1 ) 


(12) 


(13) 


DDCO 


0.805 


15.71 


1.77(43) 


1.20(43) 


2009ig 


1 


32.60 


TF 


0.14 


1.12 


-3.3 


1.03 


DDC6 


0.727 


16.89 


1.62(43) 


1.03(43) 


2004dt 


2.3 


34.44 


2 CMB 


0.04 


1.20 


-1.9 


1.08 


DDC10 


0.650 


17.44 


1.45(43) 


1.12(43) 


2005cf 


4,5 


32.17 


SNIa 


0.22 


1.00 


-0.8 


1.05 


DDC15 


0.558 


17.76 


1.22(43) 


9.25(42) 


2002bo 


6 


31.90 


SNIa 


0.41 


1.06 


-0.2 


1.05 


DDC17 


0.458 


18.90 


1.03(43) 


8.51(42) 


1992A 


7 


31.63 


SBF 


0.02 


0.98 


+2.7 


0.91 


DDC20 


0.344 


19.38 


7.65(42) 


4.15(42) 


2()07hj 


8 


33.88 


■^Virgo 


0.11 


1.32 


-2.0 


1.28 


DDC22 


0.236 


19.52 


5.10(42) 


3.39(42) 


2003gs 


9 


31.65 


SBF 


0.07 


1.12 


+1.8 


1.19 


DDC25 


0.180 


20.50 


3.64(42) 


1.91(42) 


1999by 


10 


30.74 


Cepheids 


0.03 


1.38 


+2.7 


1.47 



Notes: Numbers in parentheses correspond to powers of ten. Column information: (1) Model name; (2) mass of synthesized 56 Ni; (3) rise time to peak 
bolometric luminosity. Note that the rise to UVOIR maximum is ~0.5— 1 d longer than to t he true bolometric m aximum; (4) peak bolometric luminosity; (5) 
peak UVOIR luminosity derived from synthetic U BVRI magnitudes using the method of lValenti et al.l2008T ; (6) name of observed SN la used for 
comparison in Fig. [5] (7) reference for spectroscopic data used in Fig.|5](see References below); (8) distance modulus (assuming Ho = 73kms -1 Mpc -1 
where applicable); (9) method for distance modulus determination: SBF=Surface Brightness Fluctuations, SNIa=SN-based distance, TF=Tully Fisher, 
2CMB=CMB-frame redshift-based distance modulus; zvirgo=distance modulus based on the recession velocity correcte d for the influence of the Virgo 
cluster ; (10) total (Galactic and host-galaxy) reddening. Galactic reddening values are based on the infrared dust maps of lSchlegel. Finkbeiner & Davisl 
119981) . while host-galax y reddening was infer red based on our own fits to optical (and NIR when available) SN lightcurves. For SN 2003gs we use the 
reddening reported bv lKrisciunas et ai]d2009h : (11) ra tio of synthetic to o bserved peak UVOIR luminosities, derived from UBVRI magnitudes (only 
BVRI in the case of SN 2007hj) using the method of lValenti et al.l2008l ; (12) phase of SN with respect to the time of UVOIR maximum; (13) scale factor 
applied to the ob served spectra to matc h the mean synthetic flux in the wavelength range 4000-6000 A. 

Re ferences: (l):lFolev et al.l2012l (2):lAltavilla et al.l2007t (3): IWang et al.l2012l ( 4): lGaravini et al.l2007t (5): lBufano et al.l2009l (6): lBenetti et al.l2004 
(7): lKirshner et al.ll993l (8): ISilverman et al.l2012t (9): Kotak et al. (in prep.), (10): lGarnavich et 31.120041 



The broader S i II lines in the spectra of SN 1997bp 
jBlondin et ai1l2012h are better matched by DDC10, but this SN 
only has published spectra in the range 3500-7500 A. 

6.4 DDC15 (M 5 6 Ni = 0.56 M Q ) vs. SN 2002bo 

There are only SN-based distances to SN 2002bo, which vary by 
~ 0.5 mag in distance modulus (i.e. ~ 20% in absolute distance). 
This uncertainty stems in part from the large reddening inferred 
for this SN (E(B — V) ~ 0.4 mag, the largest of our comparison 
sample; see Table [3} . We adopt the distance modulus of 3 1.90 mag 
inferred from NIR light curves dWood-Vasev et alj|2008l) . less sub- 
ject to extinction than the visual bands. A similar distance mod- 
ulus is derived from NIR light curves of the highly-reddened 
SN 2002cv (/j = 31.85 mag assuming H = 73kms _1 Mpc -1 ; 
IWood-Vasev et al]|2008l) . which occurred in the same host galaxy 
as SN 2002bo. Our inferred host-galaxy reddening (E(B — V) = 
0.39 mag from our own fits to the lightcurves using th e BayeSN 
statistical model of iMandel. Narayan & Kirshneril201 ll) is consis- 
tent with published estimates (e.g. jBenetti et al.l2004l) . The authors 
of this latter paper argue based on spectroscopic modelling that a 
lower value of the reddening (E(B — V) ~ 0.30 mag) is needed to 
match the maximum-light spectrum, although their synthetic spec- 
tra systematically overestimate the flux redward of ~ 6500 A due 
to their assumption of an inner diffusive boundary (their Fig. 17). 

The synthetic spectrum of model DDC15 and the observations 
of SN 2002bo nearly overlap. The peak UVOIR luminosity is well 
reproduced (Q uvo ir = 1.06, i.e. better than 10%), as are the mag- 
nitudes at bolometric maximum in the optical and NIR bands (see 
Table CI). 

W e obtain equally-good matches to SN 2002dj jPignata et alj 
20081). SN 2004fu. a nd SN 2007co (both with spectra published by 
Blondin et ai]|2012h . all of which are also Broad Line SNe la of 
comparable luminosity. 



A decent mat ch is also found to the nearby SN 201 lfe 
dNugent etai]|201ll) . although our model DDC15 is significantly 
redder at maxi mum light (but of comparable luminosity, with 
Quvoir ~ 1.1). iRQpke etalJ d2012h found a reasonable match to 
SN 20 life using a 3D delayed-detonation model tuned to synthe- 
size 0.6Mq of 56 Ni, although they argue that the agreement at max- 
imum light is better for their violent WD merger model. 

DDC15 is al so a good match to SN 200 lay, for which 
iBaron et"al] d2012h invoke a pulsating delayed detonation model 
with a special WD configuration (excess C abundance and low 
central density). Our standard delayed-detonation model DDC15 
is able to reproduce the peak UVOIR luminosity (Quvoir = 1.05) 
and colours, as well as the broad spectroscopic features (in partic- 
ular Si II 6355 A). 

6.5 DDC17 (Ms6 Ni = 0.46 M ) vs. SN 1992A 

iKirshner et ai] d 19931) present an extensive set of UV spectra for 
SN 1992 A taken with the International Ultraviolet Explorer (IUE) 
as well as HST. The maximum-light spectrum shown in Fig. [5] is a 
combination of an IU E UV spectrum and a CTIO optical spectrum. 
IKirshner et al. I< fl993h obtain reasonable fits to the observed spectra 
using a parameterized steady-state LTE radi ative transfer ca lcula- 
tion of the delayed-detonation model DD4 of |WoosleyJ dl99ll) . Our 
full non-LTE time-dependent calculation reproduces the spectrum 
over the entire observed wavelength range, and is consistent with 

the inferred peak UVOIR luminosity (Quvoi r = 0.98). 

SN 2004eo (labeled "transitional" by IPastorello et alj|2007h 
is also a fair match to model DDC17, although Q UV oir ~ 0.7, 
i.e. model DDC17 is fainter. Yet its 56 Ni mass (046M S ) is ~ 
0.1 Mq larger than that inferred bv lMazzali et alj d2008l) based on 
modelling a s i ngle n ebular-epoch spectrum (0.38 Mq). However, 
iMazzali et alj d2008h note that their inferred 56 Ni mass falls short 
of reproducing the peak bolometric luminosity, and revise their es- 
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timate to 0.43 ± 0.05 M , in agreement with model DDC17. A 
clear asset of our approach is the simultaneous calculation of the 
light curves and the spectra, which ensures the overall spectral mor- 
phology accurately reflects the SN luminosity. 



6.6 DDC20 (Ms. 



0.34 M ) vs. SN2007hj 



Model DDC20 corresponds to a transition from standard to low- 
luminosity SNe la, with the noticeable impact of increasing blan- 
keting from Ti Il/Sc II around ~ 4000 A. 

A good match is found to SN 2007hj. This SN is slightly bluer 
than the model but the relative strengths of spectral features is well 
reproduced, in particular the 4000-4500 A range, the S II "W" fea- 
ture, the Si II lines (5972 A and 6355 A) and corresponding 7£(Si) 
ratio, the O 1 7773 A absorption, the Ca II near-IR triplet, as well as 
the Mg II feature at ~ 8800 A (see Fig. B6). There are however sig- 
nificant differences in the blue wing of the Ca II H&K line which 
is strongly suppressed in the model mainly due to Ca II absorption 
out to high velocities (see Sect. 15.4b . Lines of Ti II and Sc II also 

deplete the flux in this region. 

The spectroscopic features of SN 1998 bp l lMatheson et alj 

l2008h and SN 2007au l lBlondin et ail |2012|) are well matched 
by model DDC20, but the predicted luminosity is too high 
(Quvoir = 1-91 and 1. 62, respectively). The match to SN 2000dk 
l Matheson et alj |2008|) is of comparable quality as the match 
to SN 2007hj, although the discrepancy in the blue wing of 
Ca II H&K is even larger (the luminosity is however in better agree- 
ment, with Quvoir = 0.85). 



6.7 DDC22 (Af56 Ni = 0.24 M ) vs. SN 2003gs 

With DDC22 we enter the category of subluminous SNe la, with 
the characteristic broad Ti Il/Sc II absorption feature over the range 
4000^1500 A. 

An impressive match is found for SN 2003gs, a low- 
luminosity SN la in the optical bands (though slightly more lumi- 
nous than 91bg-like SNe la). We assume the time of UVOIR max- 
imum to coincide with the first photometric measurement (which 
in turn corresponds to the inferred time of l/-band maximum by 
iRrisciunas et alj|2009l) since there are no pre-maximum observa- 
tions for this SN. We expect the impact on the derived ratio of peak 
UVOIR luminosities (Quvoir = 1.12) to be minor, given the coin- 
cidence within 1 d of the times of UVOIR and V-band maxima in 
model DDC22. We checked t hat using the corr e cted t /-band pho- 
tometry recently published by iKrisciunas et al] j2012l) had no no- 
ticeable impact on the derived Quvoir- 

While the flux in the blue wing of Ca II H&K compares better 
than for model DDC20 versus SN 2007hj, the predicted HVFs in 

the Ca II near-IR triplet are not present in the obse rvations. 

The spectroscopic features of SN 1998de l lMatheson et alj 

120081) are well matched by model DDC22, but the predicted lumi- 
nosity is much too high (Q UV oir = 2.5 based on BVRI photometry 
only). 



6.8 DDC25 (Ms6 Ni = 0.18 M ) vs. SN 1999by 

Model DDC25 corresponds to the classical low-luminosity 91bg- 
like SNe la. The overall SED and individual spectral features 
match up nicely with the well-observed 91bg-like SN 1999by 
fcjarnavich et alj |2004|) . We encounter the same problem as in 
model DDC22 concerning the prediction of prominent HVFs in 



the Ca II near-IR triplet. The model appears slightly over-luminous 
with respect to this SN (Q UV oir = 1.38), but using only a la larger 
Cepheid distance to the host galaxy NGC 2841 ( i.e. fj, — 30.97 mag 
instead of — 30.74 mag; iMacri et al.l [200lh results in a better 
agreement, with Q uvo ir ~ 1.1 (i.e. 10% shift in absolute flux). 
This model compares w ell to other 91bg-like SNe la, such as the 
eponymous SN 1991 bg dFilippenko et al.lll992l ; iLeibundgut et alj 
1 1993b and SN 2005b! dTaubenberger et alj|2008h. 



We thus confirm the result of IHoflich et ail d2002h . that 
low-luminosity SNe la can be explained within the framework 
of Chandrasekhar-ma ss delayed-detonati on models, contrary to 
claims otherwise (e,g,. |Pakmor et al I20T l]). We will need, however, 
to see how our delayed-detonation models fare at pre-maximum 
phas es through to t he neb ular phase. 

iHowell etalJ feOOlh show evidence for a significant intrin- 
sic polarization level of ~0. 3-0.8% in SN 1999by near maximum 
light, and argue that this SN has a well-defined symmetry axis, pos- 
sibly resulting from a rapidly-rotating single WD or a binary WD 
merger process. Our excellent fits to the observed spectra with a ID 
radiative transfer code suggest that such asymmetries, if present, do 
not dramatically affect the observables around peak bolometric lu- 
minosity. 



7 DISCUSSION AND CONCLUSION 

In this paper, we have presented a first effort at modeling SNe la 
with CMFGEN. Our ID approach solves the radiative-transfer prob- 
lem for SN la ejecta with a high level of physical consistency, 
which comes at a great computational expense. The coupled set 
of time-dependent statistical equilibrium, gas energy, and radiative 
transfer (0 th and 1 st moments) equations are solved for simultane- 
ously, yielding multi-epoch SN la spectra from the UV to the in- 
frared. The bolometric luminosity and the photometry are obtained 
by suitable integrations of that wavelength-dependent flux. Impor- 
tantly, line blanketing is treated explicitly and consistently — the 
algorithm makes no assumption about the thermalization character 
of the opacity as typically done in approaches that do not solve for 
the level populations directly. 

We compute r adiative properti es for a sequence of delayed- 
detonation models ( lKhokhlov|[l99lh . a mechanism in which a de- 
flagration is followed, after some delay depending on some pre- 
scribed transition density, by a detonation in a Chandrasekhar- 
mass white dwarf. This model is not new, and our inves- 
after nearly two decades of study of this 



tigation co mes after nearly 
mechanism jKhokhlovl [l99lt H 



1993 



1996 



Khokhlov, Miiller & Hoflich 
Hoflich. Khokhlov & Wheeled 1 19951: IHoflich & Khokhlov! 
Hoflich et alj 120021 : iGamezo. Khokhlov & Orarj|2005h . However, 
we present results that have a higher level of physical consistency 
than previously achieved, albeit in ID. 

Our models start at 1 d after explosion and are evolved with 
CMFGEN until beyond the bolometric peak. In this paper, we 
focus on the phase of bolometric maximum of such delayed- 
detonation models. We analyze their properties and compare them 
with publicly-available SN la data. For each model in the delayed- 
detonation sequence, best characterized by a 56 Ni mass ranging 
from 0.18 to 0.81 M , we generally find a SN la that fits the 
peak bolometric flux to within ~ 10%, and the colors to within 
~0. 1 mag. The spectral features are typically reproduced with high 
fidelity, both in strength and width. 

Remarkably, with our set of delayed-detonation models, we 
reproduce the basic morphological diversity of SN la radia- 
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tive properties, from blue, luminous, yet spectroscopically-normal 
SNe la (e.g., SN2009ig), to low-luminosity 91bg-like SNe la (e.g., 
SN 1999by). Our assessment is based primarily on optical data 
where the bulk of the SN la flux emerges, but we include UV 
data from the Hubble Space Telescope and the Swift satellite when 
available. We find that the rise time to bolometric maximum in- 
creases from ~ 15 to ~21d with decreasing Ni mass. Despite 
the c omplexity of our c alculations, and the simplicity of Arnett's 
rule d Arnettll 1 9791. 1 1 982h . we support its validity to within 10 % in 
such delayed-detonation models. Hence, at its maximum, the bolo- 
metric luminosity of SNe la is close to the instantaneous rate of 
decay energy, although both quantities are directly related only at 
nebular phases (and physically equivalent in the case of full 7-ray 
trapping). 

We infer that the spectra of SNe la at maximum light primar- 
ily reflect the change in ionization balance resulting from the mag- 
nitude (and to a lesser extent, distribution) of heating from Ni 
decay. A reduction in temperature causes a reduction of the flux in 
the blue, which is further exacerbated by a growing opacity from 
recombining ions like Ti, Sc, or Fe. This effect is particularly se- 
vere in the U and B bands for models with a 56 Ni mass less than 
~0.3 Mq, and causes the transition to 91bg-like SN la spectra. 

Reproducing the line widths is more prone to error due to 
our ID approach, and the large ejecta kinetic energy in the ex- 
plosion models, which leads to broader lines than observed in the 
bulk of SNe la. This problem seems however to affect mainly 
the Si II lines, while others (e.g., Ca IlH&K, which overlap with 
Si II 3858 A in the blue, and the S II "W" feature around 5500 A) 
match the observations. Across our model set, we reproduce satis- 
factorily the widths of most features, the near-constant strength of 
Si II 6355 A, and the strong dependence of 72.(Si) on 56 Ni mass. 

Ionization effects are particularly important for Ca II lines. In 
weaker explosions, Si II 6355 A appears narrower when the Ca II 
near-IR triplet lines follow the opposite trend and appear broader 
and stronger. While the former represents an interesting diagnostic 
of the explosion energy, the latter is more sensitive to ionization. In 
the delayed-detonation model with the least amount of 56 Ni, ion- 
ization kinks give rise to high-velocity features in Ca II near-IR 
triplet: such "structures" are spherical shells of distinct Ca ioniza- 
tion, producing optical-depth variations that affect line profiles al- 
though they would produce no polarization signature. 

In our model ejecta, there is essentially no carbon left by 
the detonation, even if it is weak. This is not a shortcoming 
of the delayed-detonation mechanism, but rather a problem with 
the assumption of spherical symm etry for the explosion phase 
dGamezo. Khokhlov & Oranl l2005h - So, while ionization is key, 
abundances can also completely inhibit the formation of certain 
lines. The O 1 7773 A line is in principle a useful probe of unburnt 
oxygen, although none of our maximum-light spectra shows O I 
absorption out to velocities where we recover the initial X(0) fs 
0.5 composition of the progenitor WD. Nonetheless, our models 
reproduce the observed trend of increasing O 1 7773 A absorption 
with decreasing luminosity at maximum light, and we obtain no O 
signature in our (ID) delayed-detonation models that have a 56 Ni 
mass > 0.4 M Q . 

Even today, the inferences on the WD ejecta mass are not 
sufficiently accurate to establish whether it is the Chandrasekhar 
mass, or more, or less. At present, bot h Mph and sub-Mph mod- 
els are compati ble with observations dHoflich & Khokhlov|[l996l 
ISim et alj|20ld) . Hence, an answer to this most fundamental ques- 
tion is still lacking. To identify a good model for a SN la, fitting 
the light curve is mandatory but it is clearly not sufficient. Fitting 



the spectra is critical, and doing so simultaneously, with the same 
physics, is an asset. However, one must ensure a suitable match to 
observations at all times, not just at maximum light as done in this 
paper. In particular, explosion models are expected to best reveal 
their differences shortly after explosion (e.g. through presence or 
absence of unburnt material) or at nebular times when the IGE core 
(and possibly the innermost 56 Ni-deficient region) is revealed. 

Moreover, one is faced with the role of multi-D effects, which 
may alter the spectral appearance size ably for different viewing an- 
gles dKasen. R6pke & Wooslevl2009l) . As we demonstrate here, as- 
suming spherical symmetry may still be suitable at certain epochs 
(e.g., bolometric maximum). However, ID models miss important 
effects that are inherently multi-dimensional, such as mixing. We 
can mimic the impact of mixin g in ID by softening composition 
gradients, either artificially (e.g. jDessart et al]|2012al) . or based on 
angular av erages built from the few 3D s imula tions that exist, such 
as those o f G amezo, Khokhlov & Oranl J2005I) . Performing multi- 
D radiative-transfer simulations would seem warranted but all such 
simulations so far ha ve dramatically simplified the transfer (e.g., 
iKromer & Simll2009l) . One should question if the gain in realism 
by going to 3D is not offset by the much simplified radiative trans- 
fer compared to ID. 

An argument often made is that a given explosion model, 
when post-processed by a radiative-transfer code, sometimes fails 
to reproduce a given observation. One can then propose that this 
initial explosion model is just not perfect. However, equally good 
interpretations are that the radiative-transfer is missing some im- 
portant process, or that the atomic data is inaccurate. The advantage 
of our strategy over such targeted studies is that if a given mech- 
anism is any good at reproducing SN la observations, we should 
be able to find a number of perfect matches — as we did. Out of a 
few hundred SNe la, some must be pointing the "right" way for our 
model, or have a suitable chemical stratification, i.e., so that what 
we adopt in our ID model is a fair representation of what we see 
averaged over all viewing angles. A viable model can fail for one 
SN la, but should not fail for the few hundred SNe la for which we 
have data today. 

With the same set of delayed-detonation models, we will now 
investigate the properties at earlier times (the brightening phase 
to the peak of the light curve) and later times (transition past the 
peak of the light curve into the nebular phase). Separately, we will 
perform similar simulations for a wider range of progenitor and 
explosion mechanisms, including pure deflagrations and pulsating 
delayed detonations. 
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APPENDIX A: MODEL ATOMS 

Table A l gives the number of l evel s (both super-levels and full lev- 
els; see Iffillier & Millerll 19981 and lDessart & Hillier|[2oTol for de- 
tails) for the model atoms used in the radiative-transfer calculations 
presented in this paper. 

Oscillator strengths for CNO eleme nts were originally taken 
from lNussbaumer & Storey! d 19831 , Il984l) . These authors also pro- 
vide transition probabilities to states i n the ion conti nuum. The 
largest source of oscillator data is from lKuruczl |2009f|; its prin- 
cipal advantage over many other sources (e.g., Opacity Project) 
is that LS coupling is not assumed. More recently, non-LS os- 
cillator strengths hav e become available through the Iron Project 
dHummeret"al]|l993l) . an d work done by the atomic-data group 
at Ohio State University jNahaj 120101). Other important sour ces 
of radiative d ata for Fe include iBecker & Butlen Jl992ll 1995al lbh. 
iNahaj Jl995h . Energy levels have generally been obtained from 
National Institute of Standards and Technology. Collisional data 
is sparse, particularly for states far from the ground state. The 
principal source for collisional data am ong low ly i ng sta tes for 
a variety of spec ies is the tabulat i on by iMendozal dl 9823b; othe r 
sources include iBerrington et al. Jl985l). ILennon et all Jl985h . 

Jnskvl f 



ILennon & Burkel dl994h Ishine & LinskyTd 1974) . iTavall l ll997al lbl). 
Zhang & Pradhan jl995al : Il995bl; Il997t). Ph otoionization data is 
taken from the Opaci ty Project l Seatonlll987l) and the Iron Project 
dHummer et alj fl993l) . Unfortunately Ni and Co photoionization 
data is generally unavailable, and we have utilized crude approx- 
im ations. Charge exchange c ross-sections are from the tabulation 
bv lKingdon & Ferland fl99^) . 



APPENDIX B: CONTRIBUTION OF INDIVIDUAL IONS 
TO THE TOTAL FLUX 

In what follows we illustrate the contribution of individual ions to 
the full synthetic spectra, compared to the same observed SN la 
spectra as in Fig. [5] The upper pannel in each plot shows the full 
model (black) and the observed spectrum (red). Labels are the same 
as in Fig. 5. The lower pannel in each plot shows the contribution 
of individual ions, computed by taking the ratio of the full spec- 
trum (Ja, ah) to that excluding all bound-bound transitions of the 
corresponding ion (Fa, less)- 

We selected ions that contribute significantly to the observed 
total flux in at least one of the delayed-detonation models at bolo- 
metric maximum. The following ions showed essentially no notice- 
able contribution in the wavelength range considered here (2000- 



Table Al. Summary of the model atoms used in our radi ative-transfer cal- 
culatio ns. The source of the atomic data sets is given in iDessart & Hillieil 
(20101 1. Nf refers to the number of full levels, N s to the number of super lev- 
els, and Ntrans to the corresponding number of bound-bound transitions. 
The last column refers to the upper level for each ion treated. The total 
number of full (super) levels treated is 13 959 (2149), which corresponds 
to 629 396 bound-bound transitions. Note that nw 2 W refers to a state with 
principal quantum number n (all I states combined into a single state), and 
spin 2. Similarly, 8z X Z refers to the n = 8 state with high 1 states (usually 
I = 4 and above) combined and spin 1 . 



4 Data are available online at http://kurucz.harvard.edu 
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10000 A): C I— II, O II, Ne II, Na I, Mg I, Al II— III, Ar III, Sc III, 
Ti III, Cr III, Mn II, and Ni III. 

These plots are mostly for educational purposes, yet reveal 
key observed trends with decreasing 56 Ni mass: increase in Ca II 
strength (with double-absorption features resulting from optical- 
depth effects alone; see Sect. 5.4); increase of 7£(Si) ratio and de- 
crease of Si II 6355 A absorption velocity and width; sudden onset 
of prominent Ti Il/Sc II blanketing in the 4000-4500 A region at 
the transition to low-luminosity SNe la; increase in strength of the 
O 1 7773 A line. 



APPENDIX C: LOG OF BOLOMETRIC LUMINOSITY 
AND ABSOLUTE MAGNITUDES AT BOLOMETRIC 
MAXIMUM 

Table CI lists the bolometric luminosity and UBVRIJHK mag- 
nitudes of our models at bolometric maximum. The last line gives 
their standard deviation amongst our model series. Note the small 
scatter in the NIR bands with respect to the optical. 
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Figure Bl. Contribution of individual ions (bottom) to the full synthetic spectrum of DDCO (top, black line), compared to SN 2009ig at —3.3 d from UVOIR 
maximum (top, red line). 
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Figure B2. Contribution of individual ions (bottom) to the full synthetic spectrum of DDC6 (top, black line), compared to SN 2004dt at — 1 .9 d from UVOIR 
maximum (top, red line). 
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Figure B4. Contribution of individual ions (bottom) to the full synthetic spectrum of DDC15 (top, black line), compared to SN 20()2bo at —0.2 d from UVOIR 
maximum (top, red line). 
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Figure B5. Contribution of individual ions (bottom) to the full synthetic spectrum of DDC17 (top, black line), compared to SN 1992A at +2.7 d from UVOIR 
maximum (top, red line). 
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Figure B6. Contribution of individual ions (bottom) to the full synthetic spectrum of DDC20 (top, black line), compared to SN 2()07hj at — 2.0 d from UVOIR 
maximum (top, red line). 
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Figure B7. Contribution of individual ions (bottom) to the full synthetic spectrum of DDC22 (top, black line), compared to SN 2003gs at +1.8 d from UVOIR 
maximum (top, red line). 
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Figure B8. Contribution of individual ions (bottom) to the full synthetic spectrum of DDC25 (top, black line), compared to SN 1999by at +2.7 d from UVOIR 
maximum (top, red line). 
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Table CI. Bolometric luminosity and absolute U BVRIJHK magnitudes at bolometric maximum in our model set. 



Model 


-^bol 


Mu 


Mb 


71 /T 

My 


Mr 


7t ! 


7t ! 

Mj 


71 J 


7t J 




[erg s _1 ] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


DDCO 


1.77(43) 


-20.16 


-19.53 


-19.51 


-19.36 


-18.67 


-18.14 


-17.41 


-17.56 


DDC6 


1.62(43) 


-20.05 


-19.44 


-19.41 


-19.31 


-18.66 


-18.18 


-17.44 


-17.52 


DDC10 


1.45(43) 


-19.91 


-19.35 


-19.30 


-19.26 


-18.69 


-18.26 


-17.54 


-17.56 


DDC15 


1.22(43) 


-19.63 


-19.20 


-19.16 


-19.18 


-18.73 


-18.33 


-17.68 


-17.64 


DDC17 


1.03(43) 


-19.28 


-19.03 


-19.03 


-19.09 


-18.70 


-18.33 


-17.78 


-17.67 


DDC20 


7.65(42) 


-18.42 


-18.75 


-18.84 


-18.91 


-18.61 


-18.34 


-17.99 


-17.78 


DDC22 


5.10(42) 


-17.34 


-17.99 


-18.51 


-18.64 


-18.47 


-18.28 


-18.17 


-17.90 


DDC25 


3.64(42) 


-16.65 


-17.40 


-18.19 


-18.36 


-18.20 


-18.07 


-18.04 


-17.72 


StdDev 


5.16(42) 


1.33 


0.76 


0.46 


0.35 


0.18 


0.10 


0.29 


0.13 



Note: Numbers in parentheses correspond to powers of ten. 

Li C, Hillier D. J., Dessart L., 2012, MNRAS, 426, 1671 

Maori L. M., Stetson P. B., Bothun G. D., Freedman W. L., Gar- 

navich P. M., Jha S., Madore B. R, Richmond M. W., 2001, ApJ, 

559, 243 

Mandel K. S., Narayan G., Kirshner R. P., 201 1, ApJ, 73 1, 120 

Matheson T. et al., 2008, AJ, 135, 1598 

Mazzali P. A. et al., 2005, ApJ, 623, L37 

Mazzali P. A., Lucy L. B., 1993, A&A, 279, 447 

Mazzali P. A., Sauer D. N., Pastorello A., Benetti S., Hillebrandt 
W., 2008, MNRAS, 386, 1897 

Mendoza C, 1983, in IAU Symposium, Vol. 103, Planetary Neb- 
ulae, D. R. Flower, ed., pp. 143-172 

Milne P. A. et al., 2010, ApJ, 721, 1627 

Nahar S. N., 1995, A&A, 293, 967 

— , 2010, Norad-atomic-data 

Nomoto K., Sugimoto D., Neo S., 1976, Ap&SS, 39, L37 
Nugent P., Phillips M., Baron E., Branch D., Hauschildt P., 1995, 

ApJ, 455, L147+ 
Nugent P. E. et al., 201 1, Nature, 480, 344 
Nussbaumer H., Storey P. J., 1983, A&A, 126, 75 
— , 1984, A&AS, 56, 293 

Pakmor R., Hachinger S., Ropke F. K., Hillebrandt W., 2011, 
A&A, 528, A117+ 

Pakmor R., Kromer M., Ropke F. K., Sim S. A., Ruiter A. J., Hille- 
brandt W., 2010, Nature, 463, 61 

Pakmor R., Kromer M., Taubenberger S., Sim S. A., Ropke F. K., 
Hillebrandt W., 2012, ApJ, 747, L10 

Parrent J. T. et al., 2011, ApJ, 732, 30 

Pastorello A. et al., 2007, MNRAS, 377, 1531 

Phillips M. M., 1993, ApJ, 413, L105 

Phillips M. M. et al., 2006, AJ, 131, 2615 

Pignata G. et al., 2008, MNRAS, 388, 971 

Pinto P. A., Eastman R. G, 2000a, ApJ, 530, 744 

— , 2000b, ApJ, 530, 757 

Pskovskii I. P., 1977, Soviet Astronomy, 21, 675 

Ropke F. K. et al., 2012, ApJ, 750, L19 

Ropke F. K., Niemeyer J. C, 2007, A&A, 464, 683 

Sauer D. N. et al., 2008, MNRAS, 391, 1605 

Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 525 

Seaton M. J., 1987, Journal of Physics B Atomic Molecular 

Physics, 20, 6363 
Seitenzahl I. R., Townsley D. M., Peng E, Truran J. W., 2009, 

Atomic Data and Nuclear Data Tables, 95, 96 
Shine R. A., Linsky J. L., 1974, Sol. Phys., 39, 49 
Silverman J. M., Filippenko A. V., 2012, MNRAS, 425, 1917 



Silverman J. M. et al., 2012, MNRAS, 425, 1789 
Sim S. A., 2007, MNRAS, 375, 154 

Sim S. A., Ropke F. K., Hillebrandt W., Kromer M., Pakmor R., 

FinkM., Ruiter A. J., Seitenzahl I. R., 2010, ApJ, 714, L52 
Taubenberger S. et al., 2008, MNRAS, 385, 75 
Tayal S. S., 1997a, ApJS, 111, 459 
— , 1997b, ApJ, 481,550 
Thomas R. C. et al., 201 1, ApJ, 743, 27 

Travaglio C, Hillebrandt W., Reinecke M., Thielemann E, 2004, 

A&A, 425, 1029 
Tully R. B., Fisher J. R., 1988, Catalog of Nearby Galaxies 
Utrobin V. P., Chugai N. N., 2005, A&A, 441, 271 
Valenti S. et al., 2008, MNRAS, 383, 1485 
Wang L. et al., 2003, ApJ, 591, 1 1 10 
Wang X. et al., 2009a, ApJ, 699, L139 
— , 2009b, ApJ, 697, 380 
— , 2012, ApJ, 749, 126 
Webbink R. E, 1984, ApJ, 277, 355 
Whelan J., Iben I. J., 1973, ApJ, 186, 1007 
Wood-Vasey W. M. et al., 2008, ApJ, 689, 377 
Woosley S. E., 1991, in American Institute of Physics Conference 

Series, Vol. 232, Gamma-Ray Line Astrophysics, Durouchoux 

P., Prantzos N., eds., pp. 270-290 
Woosley S. E., Fowler W. A., Holmes J. A., Zimmerman B. A., 

1978, Atomic Data and Nuclear Data Tables, 22, 371 
Woosley S. E., Kasen D., Blinnikov S., Sorokina E., 2007, ApJ, 

662, 487 

Yaron O., Gal- Yam A., 2012, PASP, 124, 668 
Zhang H. L., Pradhan A. K, 1995a, A&A, 293, 953 
— , 1995b, Journal of Physics B Atomic Molecular Physics, 28, 
3403 

— , 1997, A&AS, 126, 373 



© 2012 RAS, MNRAS 000.mi25l 



